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Abstract. The shape of homogeneous, smooth convex bodies as described 
by the Euclidean distance from the center of gravity represents a rather re- 
stricted class Mc of Morse-Smale functions on § 2 . Here we show that even 
Mc exhibits the complexity known for general Morse-Smale functions on § 2 
by exhausting all combinatorial possibilities: every 2-colorcd quadr angulation 
of the sphere is isomorphic to a suitably represented Morse-Smale complex 
associated with a function in Mc (and vice versa). We prove our claim by an 
inductive algorithm, starting from the path graph P2 and generating convex 
bodies corresponding to quadrangulations with increasing number of vertices 
by performing each combinatorially possible vertex splitting by a convexity- 
preserving local manipulation of the surface. Since convex bodies carrying 
Morse-Smale complexes isomorphic to P2 exist, this algorithm not only proves 
our claim but also defines a hierarchical order among convex solids and general- 
izes the known classification scheme in |35] , based on the number of equilibria. 
Our expansion algorithm is essentially the dual procedure to the algorithm 
presented by Edelsbrunner et al. in |19) . producing a hierarchy of increasingly 
coarse Morse-Smale complexes. We point out applications to pebble shapes. 



1. Introduction 

In this paper we prove that the arrangement of equilibrium points on the sur- 
faces of convex solids exhausts all combinatorial possibilities and also defines a 
hierarchical order among such objects. 

The shape of homogeneous, smooth convex bodies K £ K, as described by the 
Euclidean distance from the center of gravity represents a rather restricted class 
M.C of Morse-Smale functions on § 2 . We assume that the the critical points of these 
functions are generic singularities; namely, that the Hessian of the function at such 
a point is not zero. These points are the vertices of the Morse-Smale complex (cf. 
[30]) associated with K, the edges of the complex are special integral curves of 
the gradient: heteroclinic orbits connecting singularities. In physical terms, these 
singularities are identical to static equilibrium points of K when balanced on a 
horizontal plane under uniform gravity. The study of static equilibria was initiated 
by Archimedes [25] and attracts even current interest (cf. [TT], [12], [13], [21] or 
|36|). Since the numbers of stable, unstable and saddle points (S, U, H, respectively) 
are related by the Poincare-Hopf formula as 



(1) S + U-H = 2, 
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the so-called primary classification system in [35 could be set up by omitting sad- 
dles, based on the integers {S, U} and it was shown that this system is complete, 
i.e. there are no empty classes. Our current motivation is to go beyond this result 
and establish the completeness according to a more refined (secondary) classifica- 
tion based on the topological arrangement of generic equilibria of convex bodies, 
i.e., to show that every combinatorially possible arrangement physically also exists. 

The class of general Morse-Smale functions A4 on § 2 is known to exhaust all 
combinatorial possibilities. Similarly to the primary equilibrium classification, sad- 
dles can again be omitted from the Morse-Smale complex to obtain the so-called 
quasi-dual representation [18] : every Morse-Smale complex associated with a func- 
tion in M. is isomorphic to a 2-colored (S, U) quadrangulation in class Q ([H], [19] 
and [37] ! f° r more detail see Section [2] and Figure [1]) and the inverse statement 
is also relatively easy to see by direct construction. We introduce the shorthand 
notation Q(K) for the quasi-dual representation of the Morse-Smale complex asso- 
ciated with the convex body K. Our main goal is to show that even M.c exhibits 
the same complexity as M.: 

Theorem 1. For every Q* € Q there is a homogeneous convex body K* , with a 
C°° -class boundary, such that Q(K*) = Q* . 

We introduce Q n C Q for the 2-colored quadrangulations of the sphere with S + 
U = n vertices and the class of convex bodies with exactly S+U = n extremal points 
will be denoted by JC n C JC. It is relatively easy to see (as we show in Section[2]based 
on results from [5], [5] and [35]) that by subsequent applications of an operation 
called face contraction (where two diagonal vertices in a quadrangular face are 
merged and the resulting double edges are also merged) every quadrangulation 
Q* G Qn can be collapsed onto the stem-graph P 2 (the path graph with 2 vertices 
[23| ) via a sequence of (n — 1) graphs. If we list the latter in reverse order, we 
obtain what we call combinatorial expansion sequences: 

(2) Q 2 = P2,Q3, Qi,--, Qn-i,Qn = Q* such that Qi e Qi- 

Subsequent elements of ([2]) are connected via vertex splittings (cf.[4], [9]), the inverse 
operation of face contraction. Both face contraction and vertex splitting can be 
carried out on any Q g Q n , n > 2, and thus the sequences ([2]) can be generated both 
forward (vertex splittings) and backward (face contractions). Neither the forward 
nor the backward sequence is unique; the number of combinatorial possibilities as 
well as alternative representations in other graph classes are described in Section 
[2] To prove Theorem Q] in Section [3J we show that each combinatorially possible 
vertex splitting can be realized on convex bodies by a convexity-preserving, local 
manipulation of the surface (consisting of two consecutive, local truncations) which 
we call equilibrium splitting . By this operation we can create geometric expansion 
sequences 

(3) K 2 , K 3 ,K 4 , K n ^,K n ee K* such that K % e /Q, 

which are generated in such a way that they are linked to an arbitrarily, a priori 
given combinatorial expansion sequence ^ via 

(4) Q(Ki) = Qi. 

Combinatorial sequnces belonging to an arbitrary graph Q* £ Q n can be generated 
by running ([2]) in reverse order, arriving at P 2 . Since convex bodies in class K,2 
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carrying Morse-Smale complexes homeomorphic to Pi exist (cf. }35j). subsequently 
we can run ([3]) forward to obtain the desired convex body K* with Q(K*) = Q*. 
This algorithm not only proves Theorem [T] but also defines a hierarchical order 
among convex solids and generalizes the primary scheme of |35j , which uses only the 
numbers S, U of stable and unstable equilibria to classify convex shapes. Recently, 
in [2 8) it was shown that an inductive algorithm based on |35] is incapable of 
generating all combinatorial possibilities. We also remark that unlike combinatorial 
sequences @, our geometrical expansion sequences ([3]) can be run only forward 
from an arbitrary initial convex body K. Nevertheless, in a different setting, the 
geometric analogy of face contraction was considered in [8] . 

Our expansion algorithm is essentially the dual procedure to the one presented 
by Edelsbrunner et al. [15] . producing a hierarchy of increasingly coarse Morse- 
Smale complexes (i.e. running combinatorial expansion sequences in reverse order 
compared to f2|)): a topic attracting current interest in computational topology (cf. 
[5] or [24] ). Beyond illustrating that Mc exhausts all possible combinatorial pos- 
sibilities and offering a modest link between Morse theory and convex geometry, 
geometric expansion sequences defined in equation ([3]) appear to be the natural 
building blocks for the description of pebble shape evolution under collisional abra- 
sion. 

We describe graph expansions based on earlier results in Section [3J the main 
proof of our statement is contained in Section [3J where we show the existence of 
equilibrium splittings associated with arbitrary vertex splittings. We illustrate our 
results and discuss some related issues (including pebble abrasion) in Section [4] 

2. Graph expansions 

In this section, based on earlier results, we summarize the combinatorial prop- 
erties of Morse-Smale complexes associated with general Morse functions in M. . In 
Subsection 12.11 we describe three alternative graph representations, one of which 
is 2-colored quadrangulations in class Q. Subsection 12.31 describes the inductive 
algorithm which generates Q from a single stem-graph. By isomorphism we mean 
embedding-preserving isomorphism (cf. TO] ) and by quadrangulation we mean a 
finite graph embedded on the 2-sphere such that it may have multiple edges and 
each face is bounded by a closed walk of length 4 (cf. [T] or [5]). Also, we follow 
Archdeacon et al. pQ and regard the path graph (cf. 23 J P3 with 3 vertices as a 
quadrangulation. For clarity we now list the four different classes of quadrangula- 
tions appearing in the paper: 

• <2 0,s is the class of simple quadrangulations without coloring. 

• Q° is the class of quadrangulations without coloring. 

• Q is the 2-colored version of the latter. Note that every graph in Q° is 
bipartite; that is, 2-colorable (cf. [I] and [31]). 

• Q 3 * is the class of special 3-colored (S, U, H) quadrangulations with deg(v) = 
4 for any v £ H and S + U - H = 2 (cf. Q). 

In addition, we also use the class T 3 * of special 3-colored (S, U, H) triangulations 
with deg(v) — 4 for any v € H. In this case ([1} follows from the previous properties. 
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2.1. Graph representations of Morse- Smale complexes in M. Dong et al. 
[18] introduced three different ways for representing a Morse-Smale complex associ- 
ated with a function in A4 . It was shown in (TH] and [37] that such a Morse-Smale 
complex can be uniquely represented by a 3-colored quadrangulation in Q 3 *, where 
colors correspond to the 3 types of non-degenerate critical points (maxima, minima 
and saddles), each quadrangle is bounded by a closed walk consisting of a max- 
imum, a saddle, a minimum and a saddle in this order around the face and the 
degree of every saddle is 4. Wc follow Dong et al. [IB] and call this representation 
of the Morse-Smale complex the primal Morse-Smale graph. These graphs are spe- 
cial 3-colored quadrangulations, since vertices corresponding to one color (saddles) 
are always of degree 4. These special vertices can be removed from the primal 
Morse-Smale graph without losing information (cf. |18j): first we connect maxima 
and minima in the quadrangles, resulting in a triangulated Morse-Smale graph in 
class T 3 *, then we cancel saddle points and edges incident to them (see Figured]). 
Since the degree of every saddle is 4, the resulting graph is a 2-colored quadrangu- 
lation in class Q . We call this representation the quasi- dual Morse-Smale graph 
(cf. |18j). In this section we rely on the latter, however, in Section [3] the primal 
graph representation is preferred. We remark that all three representations (pri- 
mal, triangulated and quasi-dual) are equivalent in the sense that they are mutually 
uniquely identified. 
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(5) maximum 
(unstable) 
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integral curves connecting a maximum/minimum and a saddle 
integral curves connecting a maximum and a minimum 



Figure 1. Different representations of a Morse-Smale complex, 
a) Primal Morse-Smale graph in class Q 3 * b) Triangulated Morse- 
Smale graph in class T 3 * c) Quasi-dual Morse-Smale graph in class 

Q 



Due to a well-known bijection |21] between planar maps and 2-colored quad- 
rangulations of the sphere, Morse-Smale complexes can be also represented by the 
even more compact class of planar maps as well. (A planar map is the equivalence 
class of topologically equivalent drawings of a planar graph on the sphere). The 
transformation between the two classes is shown on Figure [2] starting with a graph 
in class Q, first we connect with an edge the two (not necessarily distinct) stable 
vertices of a face, then we cancel unstable points and edges incident to them. The 
resulting graph has S vertices, U faces and H edges, thus, in this representation, the 
Poincare-Hopf formula ([T]) is equivalent with Euler's formula. Although compact- 
ness is a great advantage of this planar map representation, it is also a disadvantage: 
connection between combinatorial and geometric expansions would be difficult to 
interpret, therefore we rely on the more expressive quasi-dual representation in this 
section. 




Figure 2. Bijection between planar maps and 2-colored quadran- 
gulations of the 2-sphere. a) 2-colored quadrangulation in class 
Q b) connecting stable points on each face c) cancelling unstable 
points results in a planar map 



2.2. Face contractions and vertex splittings on simple, non-colored graphs 

in Q°' s . Vertex splitting can be defined as the inverse of a combinatorial operation 
called face contraction (cf. pQ, [9] or [32]). Let us regard Q^ s e where n 

denotes the number of vertices and let F = (vi, 1)2, V3, V4) be a face of Q^ s with 
boundary walk v\e\V2C2V-i^'i'V\G\v^^ where V\ = v$ and edges are = {Uj, Uj+i}, 
i = 1, 2, 3, 4. The contraction of F at vertices {v\, V3} produces a new quadrangu- 
lation Qn—\ € Q n 'f.i in the following way: it shrinks the region of F by identifying 
the vertices v\ and V3, the edges ei and 62, and the edges e3 and e±. Since the faces 
of Qn-x are identical to the faces of Q°' s except F, Q^-i is also a quadrangulation, 
if n > 3. (cf. Figure [3] a)). The exceptional case when Q®' s — P3 and ' s = P2 
will also be treated as a special face contraction, when e± = e&, e2 = and V2 = V4 
(cf. Figure [3] d)), although P 2 is clearly not a quadrangulation. Bagatelj [4], and 
also Negami and Nakamoto [35] defined face contraction inside the class Q 0,s of 
simple quadrangulations without coloring and showed that Q° s can be generated 
inductively from P2 via a combinatorial expansion sequence of type ([2]), consisting 
of vertex splittings, the inverse transformation of face contraction. 

2.3. Multiple edges, multiple colors: inductive generation of Q by vertex 
splitting. The quasi-dual representations Q G Q of Morse-Smale complexes asso- 
ciated with convex bodies are two-colored (S, U) quadrangulations and might have 
multiple edges, and thus a slight generalization of Subsection 12.21 is needed. On 
multigraphs Q° e Q° the named vertices and edges may not be distinct, however, 
the definition of quasi-dual representation (i.e. that these graphs have been gener- 
ated by removing the saddle points of degree 4 from the triangular representation) 
admits only two kinds of coincidences: two diagonal vertices i>2 and W4 may coincide, 
and in this case the edges e2 and ej, may coincide. These two cases are illustrated 
in Figure [3] b) and c) and we can see that in these cases both face contractions 
(and their inverse, vertex splittings) can be applied in complete accordance with 
the original definition on simple graphs in Subsection [221 Note that in Figure [3] b) 
the internal domain bordered by the edges e2 and is not a quadrangular face but 
contains additional vertices. Extension of the expansion scheme to colored multi- 
graphs is similarly straightforward. On each face diagonal vertices have identical 
color, so after generating the color-less graph Q° € Q° the 2-colored version Q e Q 
can be generated by applying two alternative coloring schemes, so if Q° can be 
generated from P 2 by vertex splittings, so can be Q. Thus we showed that the 
combinatorial expansion sequences ([2]) described in the Introduction are sufficient 
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Figure 3. Face contraction on a graph in class Q. Like in [9], 
triangles incident to some vertices indicate that one or more edges 
may occur at that position around the vertex. Here, Vi and V3 are 
minimum points; their counterparts, which remove a maximum 
point from the graph, can be performed by switching the colors. 



to generate the quasi-dual representation of any Morse function associated with a 
convex body K £ JC. Conversely, given any such a body, the associated sequences 
can be generated by running ([2]) backward via face contractions. 

Since a face contraction can be applied to any face of a quadrangulation, the 
sequence of graphs from Q n £ Q n to Q2 — P2 is, in general, not unique. In special 
applications additional criteria may be applied to single out one sequence among 
the combinatorially possible ones. Face contraction on multigraphs was already 
used in [19] to simplify Morse-Smale complexes. Edelsbrunner et al. used the 
primal representation of the complex in which a face contraction (defined in the 
quasi-dual representation) emerges as a double edge contraction |14j . Their goal 
was to produce a hierarchy of increasingly coarse Morse-Smale complexes, therefore 
they applied an extra metric criteria (growing height differences on the edges of the 
graph) which resulted in an unambiguous sequence of graphs. This algorithm was 
also used in [16] to identify macroscopically perceptible static equilibrium points of 
3D scanned pebbles. 

As here we do not use any metric criteria while applying face contractions on a 
quadrangulation, the combinatorial expansion sequence @ is not unique. Thus, in 
Section [3] we have to show that all possible vertex splittings, leading from a graph 
Qk to any Qk+i can be geometrically performed on a convex body as an equilibrium 
splitting in a way that convexity is not affected. At an arbitrary vertex Vi of Qk, 
vertex splitting increases the degree of two incident vertices of Vi (on Figure El 
the split vertex is v\ — v$, the incident vertices are V2 and V4). If t!j has degree 

d(vi), we have Pi = f ^ possibilities to choose the two (not necessarily 

distinct) neighbors. The possibilities for splitting a vertex with degree 4 are shown 
on Figured] a). Since Qk has k vertices, we have altogether Pi possibilities 
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to generate a graph Qk+i from Qk- In Section [3] we prove Theorem [T] by showing 
that all these possibilities can be geometrically constructed; that is, a convex body 
can be modified in the vicinity of an arbitrary stable or unstable point in such a 
way that it remains convex and its Morse-Smale complex is modified according to 
an arbitrary splitting of that vertex. 




Figure 4. Vertex splittings at vertex v. Here, v is a minimum 
point (of course, v could also be a maximum), a) Vertex split- 
ting in the quasi-dual representation in class Q (d(v) — 4) b) 
Geometric interpretation of the previous vertex splittings as split- 
ting equilibrium point p in the triangulated representation in class 
T 3 *(d(p) = 8). The truncated piece of the body is illustrated by 
the white domain. 
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3. Geometric realization of graph expansions 

Throughout this section, by the center of a convex body K we mean its center 
of gravity. Furthermore, the Morse-Smale complex of K is meant to be the Morse- 
Smale complex defined on bd K by the Euclidean distance function from the center 
of K . If we measure distance from a different point w, then we write about the 
Morse-Smale complex of K with respect to w. To represent these complexes, in this 
section we use the primal representation, introduced in the previous section. 

Our algorithm is based on the repeated applications of two truncations T\ and 
T 2 which are described in Subsections 13.21 and 13.31 respectively. Inside these steps, 
we also use a smoothing subroutine S which is described in Subsection l3.ll Starting 
from the mono-monostatic body Ki £ K-2 (associated with the path-graph P2 = Q2, 
cf. [35]) these steps are applied in 2-cycles T\,T% to obtain the element K i+ i £ K4+1 
from Ki £ Id in the geometric expansion series ([3]), with prescribed Morse-Smale 
representation Qi+\ £ Qj+i. 

Throughout these geometric steps, we assume that: 
K is a convex body with a C°° -class boundary and nondegenerate critical points 
such that bd K has nonzero principal curvatures at these points. 
For simplicity, we call this property (*). Recall that a critical point p is nondegen- 
erate, if bdK is of class C 2 at p, and the Hessian of the Euclidean distance function 
at p £ bd K is not zero. 

Below we briefly introduce the steps T\, T2 and subroutine S before giving a 
detailed description. 

Subroutine S: Smoothing: Here we start with a convex body K satisfying 
(*) except that it has a C 1 -class boundary and construct a C°°-class small 
perturbation K' of K, such that K' satisfies (*), with a Morse-Smale com- 
plex homeomorphic to that of K. This perturbation is carried out in such 
a way that if bd K is a spherical cap in a neighborhood of a critical point, 
then the same holds for a neighborhood of the corresponding critical point 
of K'. 

Step T\: Truncation with a 'sphere': In this step, for any K satisfying (*), 
we construct a convex body K' satisfying (*), with a Morse-Smale complex 
homeomorphic to that of K, such that the critical point of K' corresponding 
to p has a spherical cap neighborhood in bdK' . 

Step T2: Vertex splitting: In this step, for any convex body K satisfying 
(*), and any stable or unstable point p of K that has a spherical cap 
neighborhood in bd K , and also any combinatorially possible splitting of p, 
we construct a convex body K' satisfying (*), such that its Morse-Smale 
complex is homeomorphic to the 'split' complex of K, using the chosen 
splitting. 

Throughout this section, we use the following notations and assumptions. For a 
point q £ M 3 and r > 0, the open (respectively closed) ball with center q and radius 
r is denoted by B r (g) (respectively, B r (g)). For simplicity, we set B = Bi(o) and 
B = Bi(o). We consider a convex body K C M 3 with a stable or unstable point 
p £ bd K . Without loss of generality, we may assume that the center of K is the 
origin, and p = (0,0, 1). We denote the integral curves of K that connect p to a 
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saddle point of K, in cyclic order around p, by Tx, T2, ■ ■ ■ , T^ = IV Furthermore, 
we let U be the unit tangent vector of Ti at p. 

3.1. Smoothing subroutine S. Here we show that if K satisfies (*) except that 
it has a C^-class boundary, then it can be modified to a convex body K' satisfying 
(*), with a Morse- Smale complex homeomorphic to that of K. 

Observe that, as o is the center of K, we have o £ int^T. Let / denote the 
distance function of K\ that is, / is defined by f(u) = min{A : A > 0,u £ XK} (cf. 
[6]). Let h : M 3 —> R be a nonnegative, C°°-class function such that suppft. C B 
and L 3 h(x) dx = 1. Such a function is called by various names in the literature: 
mollifier (cf. |20j). or bump junction (cf. |27j ) or probability distribution. Clearly, 
we may choose h in a way that its symmetry group is O 3 . Observe that by setting 
ht(x) — *' , we obtain a family of C°°-class functions with supp/it C iB and 
/ M3 ht(x) dx= 1. 

Let us define the function : R 3 — >• M as the convolution 



Clearly, for every t, gt is C°°-class, and as the integral average of convex functions 
is convex, it is convex. In particular, it follows that the set K t = ([0,1]) is 
compact and convex, and hence it is a convex body for sufficiently small values of 
t. Note that bdK t is a C°°-class submanifold of R™. 

Furthermore, it is known (cf. [22] or [37]), that if U is a compact set and / 
is C"-class on U for some 2 < n < 00, then on U gt converges uniformly to /, 
together with its derivatives up to order n, as t — > 0. Thus, for sufficiently small 
values of t, that is, for any t G [0, e] for some e > 0, the Morse- Smale complex of K t 
with respect to o is homeomorphic to that of K] and bd K t has nonzero principal 
curvatures at the critical points of K. This implies that K t satisfies (*). 

We are left to show that, for small values of t, the Morse-Smale complex of K t 
with respect to its center is homeomorphic to that with respect to o. To do this, 
we observe that for any t <E [0,e] there is some p > such that the Morse-Smale 
complex of K t with respect to u is homeomorphic to that with respect to o, for 
any u e intB p . Let the radius of the largest such ball be pit). Clearly, p(t) is 
a continuous function of t, and thus, by compactness, there is a value of p such 
that for any t £ [0, e], the Morse-Smale complex of K t with respect to any point of 
int Bp is homeomorphic to that with respect to o. Note that the center ot of K t 
tends to o as t — > 0. Thus, if t is sufficiently small, \o t \ < p. Choosing t with this 
property and setting K 1 = K t , the assertion readily follows. 

Remark 1. Note that if the symmetry group of / and his O 3 , then the same holds 
for g t for every value of t. Thus, if p £ bdK has a spherical cap neighborhood, 
then so does the corresponding point of K' . 

3.2. Step T%: Truncation with a sphere. We show that if K satisfies (*), and 
has a stable or unstable point at p, then there is a local modification K' of K 
such that K' satisfies (*), the Morse-Smale complex of K' is homeomorphic to that 
of K, and the stable/unstable point of K' corresponding to p has a spherical cap 
neighborhood in bdK'. 
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Recall that we may assume that p = (0,0, 1) and that the center of K is the 
origin. Then the condition that p is nondegenerate is equivalent to saying that 
||p||ki ^ 1 ^ ||p||k2, where K\ and K2 are the two principal curvatures of bdK at 
P- 

Now we define a two-parameter family of truncations of K, denoted by Kr(e). 
Let z = f(x, y) be a function the graph of which is a neighborhood of p in bdK. 
Furthermore, for any <fi £ [0,7r], let = f (t cos (j>,t sin cj>) be the restriction of / 

to the line L containing (cos 0, sin </>) and the origin. Consider a value of R such 
that is strictly smaller than any of the two principal curvatures of bd K at p. 
Then < |/^'(0)| for any line L. Since bdK is C°°-class in a neighborhood of p, 
so is / in a neighborhood of (0, 0). Thus, for some 8 > 0, < f'l(t) for any \t\ < 8. 

We choose 8 satisfying this condition, and consider the part S, with points for 
the coordinates of which x 2 + y 2 < 8 and z > 0, of the sphere of radius R that 
touches the plane {z = 1} at p from the side containing the origin. Clearly, S is 
a closed spherical cap. Now, set Cr as the intersection of all closed half spaces 
that contain 5* and with their boundaries tangent to S. Then Cr is a 'cone with a 
rounded apex'. Note that Cr is convex, contains K, and that (bdCfl) D K = {p}. 
Furthermore, the boundary of the translate Cr(e), with a sufficiently small vector 
(0, 0, — e) where e > 0, intersects bd K in a simple smooth closed curve contained 
in the translate of the interior of the spherical cap S. Our first lemma describes 
the properties of the intersection Kr(e) = Cr(e) D K. 

Lemma 1. Let be smaller than any of the two principal curvatures of bdK at 
p. Then there are constants Ai,A2, A3 > depending only on R and K such that 
for every sufficiently small e > the following hold. 

(1) no point of (bdK) n B Al ^(p) belongs to Kr(e). 

(2) (bdK)\B X2 ^(p)CK R (e). 

(3) If D is a convex body, with center d and satisfying Kr(e) CDC K, then 
\o'\ < X 3 e 2 . 



Proof. First, observe that Kr(e) is convex, as it is the intersection of convex bodies. 

Recall the two- variable function / from the definition of Kr(e). Like there, let 
/\l ■ K —> K denote its restriction /|z,(i) = fit cos <f>, t sin <fi), where <j) € [0, 7r], 
on a line L passing through the origin and (cos</>, sin^). Let h denote the one- 
variable function, defined on [—R,R], the graph of which is a semicircle, of radius 
R, and with maximum h(0) = 1 — e. By the conditions in Lemma [I] we have 
> f\'[(0) > h"(0). Hence, from the second-degree Taylor polynomials of /|l and 
h, we obtain that 

e - X\t 2 < f L {t) - h(t) <e- X* 2 t 2 

for some positive constants A^A^. Thus, (1) and (2) clearly follow with Ai = jr 
and A 2 = -nr. 

Now we show (3). Let d = (x' , y' , z'). Recall that 
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First, we estimate x' . Note that the part of K outside D can be covered by 
an axis-parallel brick of side- lengths 2\2yfe, 2\i\fe and ce > for some constant c 
independent of s. Thus, 



1 dV - / 1 d^ 

qED Jq£K 



and 



< 4A 2 .c£ 2 , 
< 4X 3 2 ce 5/2 , 



xdV- xdV 

IqeD Jq£K 

from which \x'\ < A^e 2 follows. In the same way, we may obtain similar bounds for 
y' and z' , which readily yields the assertion. □ 

We note that Lemma Q] can be applied for the degenerate case R = oo as well. 
The next lemma guaranties that the truncated body Kr{e) has the same numbers 
and types of static equilibrium points. 

Lemma 2. Let p be a stable point of K , or let p be an unstable point and R < \ \p\\. 
If R satisfies the conditions in Lemma\]\ then there is an e > such that for every 
t £ [0, e] the following holds: Kji(t) and K — Kr(0) have the same number of 
stable/unstable and saddle points. Furthermore, the coordinates of these points are 
continuous functions oft. 

Proof. Note that as the Euclidean distance function from a fixed point is C°°- 
class at every point but the origin, all the partials of any order at any point of 
bd-R' change continuously when translating K. Hence, applying the Poincare- 
Hopf Theorem |2, to a compact neighborhood of any equilibrium point shows that 
the numbers of the stable/unstable/saddle points of K do not change under a 
translation by a small vector, and their coordinates are continuous functions of the 
translation vector. Similarly, the gradient vector field changes continuously under 
a translation of K . Thus, o has a neighborhood U such that for any u £ U, the 
Morsc-Smalc complex of K is homeomorphic to that of K + u, or, in other words, 
the Morse-Smale complex of K with respect to u is homeomorphic to that with 
respect to o. 

If t is sufficiently small, then the center ot of Kn(t) is contained in U . Hence, all 
the stable/unstable and saddle points of K but p change continuously as functions 
of t. Note that by (1) and (3) of Lemma [l] the stable/unstable point of K, with 
respect to Ot, that corresponds to p is contained in the part truncated by Cn(t). 
Thus, this point does not belong to Kr(£). On the other hand, by (2) of the same 
lemma, there is exactly one equilibrium point of Kn(t) on the part belonging to 
C R (t). 

We leave it to the reader to show that no plane, perpendicular to [u, q] for some 
u £ V, supports Kn(t) at a point q € (hdK) n bdCu(t), and thus, there are no 
more equilibrium points of K^(t). □ 

Now we show that Step T± can be carried out using a sphere of radius R for some 
value of R. To do this, consider some Kr(e), where e > and R > satisfy the 
conditions of Lemma [2] By this lemma, we may assume that Kr(e) has the same 
numbers of stable/unstable and saddle points as K does. Let or denote the center 
of Kr(s), and pr be the critical point of Kr(e) corresponding to p. Consider a 
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closed C 1 -class curve g C bd Kr(e), separating pr from all the other critical points 
of Kr(e), which is disjoint from the truncated, spherical part, and is not tangent 
to any integral curve of bdK , with respect to or, that intersects g. 

Now we show that even though bd Kr (e) is not C 1 -class, we may still apply 
the smoothing algorithm described in Subroutine S. More specifically, we show 
that if we smooth Kr{e) using the function h t described in the subrutine, then the 
Morse-Smale complex of the resulting body K' is homeomorphic to that of K . Let 
g' C bdK' be the curve corresponding to g. Observe that as the gradient vector 
field changes continuously as a function of t, for sufficiently small t g' is not tangent 
to any integral curve of K' intersecting it. Furthermore, any such integral curve 
ends at a critical point in the region inside g' and for small values of e and t there is 
a unique critical point in this region, which we denote by p' . Thus, we may choose 
values of e and t such that the Morse-Smale complex of K' is homeomorphic to 
that of K. By Remark [2 p' has a spherical neighborhood, and the radius of this 
sphere is arbitrarily close to R. 

As a result of our consideration, we may assume that the examined critical point 
p has a spherical cap neighborhood in bd K of radius R for some value of R. In the 
next lemma, we show that this truncation can be carried out using a sphere of any 
radius r, where < r < \ \p\ | if p is an unstable point, and r > \\p\\ if p is a stable 
point. To formulate it, recall that the integral curves connecting p to a saddle point 
of K are denoted by Ti, T2, . ■ ■ , I 1 /. = Tq in cyclic order around p. 

Lemma 3. Let K satisfy (*), and assume that p has a spherical cap neighborhood 
in bdK, of radius R. Let < r < \\p\\ if p is an unstable point, and r > ||p|| if 
p is a stable point, and let 6 > arbitrary. Then there is a convex body K' C K 
satisfying (*) and the following. 

• The Morse-Smale complex of K' is homeomorphic to that of K. 

• Denoting the critical point of K' corresponding to p by p' , p' has a spherical 
cap neighborhood in bdK' , of radius arbitrarily close to r. 

• Denoting the integral curve of K' corresponding to Ti by T^ for every i, and 
by t'i the unit tangent vector o/Tj at p' , we have that — ti\ \ < 8. 

Proof. First, let r < R. Then we set p' — p. Let U be a (sufficiently small) 
neighborhood of o such that for every u £ U, the Morse-Smale complex of K with 
respect to u is homeomorphic to that with respect to o. By (3) of Lemma [TJ there 
is an e > such that for any convex body K', satisfying (K \ B e (p)) C K' C K, 
the center of K' is contained in U. 

Clearly, we may replace the part of bd K in B e (p) with a surface of rotation S', 
with the line containing [o,p] as its axis of symmetry, in such a way that: 

• p £ S' C K and (bd K \ B £ (p)) U S' is the boundary of a convex body K', 

• the modified body K' has a C 1 -class boundary, 

• a neighborhood of p in S' belongs to a sphere of radius r, 

• there is no equilibrium point on S' but p, with respect to any point of U 
on the line containing [o,p\. 

Note that the center of K' lies on the line containing [o,p], and thus all the 
integral curves emanating from p in B e (p) are the meridians of S' . Let these curves 
be denoted by r' l5 . . . , r' fc = T' in cyclic order, with t\ denoting the unit tangent 
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vector of 1^ at p. Thus, since the integral curves of K with respect to w € U, 
connecting p and the saddle points of K are deformed continuously when modifying 
u, the assertion readily follows for sufficiently small U . To obtain a convex body 
with a C°°-class boundary, we may finally apply Subroutine S. 

If r > R, we may use a truncation by a sphere of radius r, like in the previous 
part of Step Ti. □ 

Remark 2. Observe that the two critical points of the convex body K with Morse- 
Smale complex P2 are the 'poles' of the parametrization (cf. |35j). and thus, may 
not have C°°-class neighborhoods. Nevertheless, the curves PnbdK, where P is 
any closed half plane containing a critical point of K and its center in its relative 
boundary, are C°°-class curves with nonzero second derivatives at the critical points. 
Thus, a straightforward modification of our arguments can be applied in this case. 

3.3. Step T2: Vertex splittings. In this subsection, given a convex body K sat- 
isfying (*), with a critical point p that has a spherical cap neighborhood in bd K, 
we describe the equilibrium splitting corresponding to any given combinatorially 
possible splitting defined in Section [5] such that the resulting convex body K' satis- 
fies (*). An example for a complete list of combinatorially possible vertex splittings 
is illustrated in Figure [4] a) , the corresponding approximate geometric truncations 
are shown in Figure |4]b). First, we consider the case that p is a stable point of K. 
The truncation is illustrated in Figure [5j 




Figure 5. Truncation with a plane 



Recall that the center of K is the origin, and p = (0,0,1). Furthermore, the 
integral curves of K, connecting p and the saddle points ot K, in cyclic order 
around p, are Ti, where i = 1, 2, . . . , k, and the unit tangent vector of Ti at p is U. 
Let R and q denote the radius and the center of the sphere, containing a spherical 
cap neighborhood of p. We denote this neighborhood by S. Then, according to 
Lemma[3]in Step Ti, we may assume that R is arbitrarily close to any given number 
greater than ||p||. Thus, the arc of any integral curve of K within S, including the 
r^s, are great circle arcs. 

Let G\ and Gi be two open great circle arcs in S, starting at p, that are, in 
the interiors of the two regions, bounded by the I\-s, that we want to split. These 
different regions are shaded in Figure Q] b) . We may choose G\ and G2 in a way 
that they do not belong to the same great circle. We intend to truncate K near 
p with a plane in such a way that there is a new stable point on the plane, a new 
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saddle point on its boundary, and the saddle point is connected to the stable point 
corresponding to p, to the stable point on the plane and to the two unstable points 
in the boundaries of the regions containing G\ and G2. To do this, we use a planar 
section that touches G\ and Gi. To examine the properties of these sections, we 
first prove a technical lemma that we are going to use in the construction. 

Lemma 4. Let C be a circle in the plane M 2 with the origin as its center and with 
radius 1. Let p = (0,i) with t>0. Let q\ = and qi = (1*2,^2) be two points 

of C such that vi > 0. 

(1) If [91,92] is perpendicular to [p,qi], then lim ^ = jzti- 

(2) If the angle of [91, 92] and [p, q{\ is ^ — cu\ for some constant c independent 
ofm, then lim ^ = ^ + 2c. 

Proof. We note that U2 — — ^7* = , from which (1) immediately follows. The 

second part can be proven with a similar elementary computation. Note that in 
the second case the angle between [p, q{\ and the segment connecting p and its 
orthogonal projection on [91,92] is cu\. □ 

Now, consider a spherical circle on S that touches G\ and G2, and let u\ and 
U2 denote, respectively, the distance of its closest and its farthest point from the 
z-axis. Let L denote the limit of 1*2/1*1 as u\ tends to zero. Note that this limit 
exists, is positive and does not depend on the radius of S, only on the angle of G\ 
and G2; this follows from the observation that spherical space 'locally' is Euclidean. 

By Lemma [3J we may slightly deform K in a way that G\ and G2 are still in 
the interiors of the two chosen regions, and that 2(i2 — 1) < L. On the other hand, 
by (1) of Lemma[4j if 2(R — 1) < L, then there is a sufficiently small circle G on S 
touching Gi and G2 such that the truncation of K by the plane containing C has a 
stable point, with respect to o, on the truncated, disk-shaped part. Furthermore, by 
(2) of the same Lemma, we may choose < r < 1 in a way that if C is sufficiently 
small, then the distance of this point from the boundary of C is at least r times 
the radius of the circle. 

Let H E denote the plane perpendicular to and intersecting [o, p] , at the distance 
e from p. Let K e denote the truncation of K by this plane. Let U e be the set 
of the centers of all the convex bodies D satisfying Ki e C D C K. Clearly, we 
may choose a sufficiently small e such that the Morse-Smale complex of K with 
respect to any point of U e is homeomorphic to that with respect to o. Furthermore, 
using (3) of Lemma [U we may choose e in a way that both the diameter of U £ , 
and the diameter of the projection on S of U E from 9 is at most ce 2 for some 
c > independent of e. Note that for small values of e, these projections belong 
to S. Furthermore, for any fixed <5 > 0, we may also choose e in a way that the 
intersection points, with the circle K n H £ , of any integral curve connecting the 
projection of a point u from 9 on S and a saddle point of K with respect to u 
moves on an arc of angle not greater than 5, where u runs over U e . We denote 
these integral curves by T%, T%, . . . , T% = T%. 

Let C £ be the circle that touches G\ , G2 and H e on the side of H £ not containing 
o, and K e be the truncation of K by the plane containing C £ . Note that the radius 
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of C £ is at least cy/e for some positive constant c. Thus, for sufficiently small e, the 
following are satisfied. 

• for any u G U e , C £ intersects exactly those curves from amongst the Tfs 
that we want to intersect with the planar section. 

• No such intersection point is on G\ and G2. 

• The orthogonal projection of U e on the plane of C e is contained in the 
interior of C e . 

• Ug is contained in K e , and is disjoint from C e . 

Clearly, under the conditions described in the previous paragraph, K E has three 
equilibrium points in the closed half plane bounded by H e and containing p: a 
stable point p 1 'near' p (that is, in U^), another stable point s on the planar disk 
bounded by C E , and a saddle point h on C e . This saddle point is connected to both 
p' and s: the integral curve connecting it to p' is a great circle arc, and the other 
one is a straight line segment. Furthermore, for every T,; on the truncated side of 
Gi U G2, there is a (piecewise differentiable) integral curve of K e connecting s to 
the corresponding saddle point of K e , and for every I\j on the other side of G\ U G2, 
there is a similar curve ending at p' . This implies, by exclusion, that h is connected 
to the two unstable points in the two chosen regions containing G\ and Gi ■ 

Finally, observe that by replacing the truncating plane containing C e by a ball 
of sufficiently large radius results in a homeomorphic Morse-Smale complex. Fur- 
thermore, 'to smooth out' bdK E , we may replace C e by a sufficiently small toroidal 
surface, which results in a convex body K' with a C -class boundary, such that its 
Morse-Smale complex is homeomorphic to the one splitted in the required way. We 
remark that Ki e C K', and thus, the center of K' is still contained in U £ . Finally, 
we may apply Subroutine S to obtain a convex body with a C°°-class boundary, 
which finishes the proof for a stable point. 

For an unstable point we may apply a similar consideration. In this case, instead 
of a truncation with a ball of large radius, we expand the original ball with a circular 
cone the axis of which contains the ball center q. 

4. Summary and related issues 

The central idea of our paper is to associate vertex splittings with localized geo- 
metric transformations. The latter are defined in such a way that we can control 
all combinatorial possibilities. Next we show that vertex splittings arise in a spon- 
taneous way in various geometric settings where they may or may not exhaust the 
full combinatorial catalogue. From this point of view, our construction creates a 
framework to study these processes. 

4.1. A roadmap for the generic bifurcations of one-parameter vector 
fields. Generic 1-parameter families of vector fields on the 2-sphere produce two 
types of singularities: saddle-node bifurcation (which results in a vertex splitting 
or a face contraction on the Morse-Smale complex) and saddle-saddle connection 
[3j- Each vector field can be uniquely associated with the quasi-dual graph repre- 
sentation of its Morse-Smale complex, so the evolution of one-parameter families 
can be studied on a metagraph Q the vertices of which are graphs Q £ Q, rep- 
resenting the Morse-Smale complexes and the edges of Q correspond to generic 
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bifurcations in one-parameter families. Any such family will appear as a path on 
Q. A small portion of Q is illustrated in Figure |6)/al. Vertices are classified based 
on the number of stable («S) and unstable (U) equilibrium points, solid edges rep- 
resent vertex splittings, dashed edges represent saddle-saddle connections, Figure 
[6]/a2 shows the latter inside the primary equilibrium class {S,U} = {2, 3}. Convex 
bodies associated with some selected graphs (selected vertices of Q) are illustrated 
in Figure |6^a3. Observe that Q is not oriented, new critical points may emerge 
or disappear at generic bifurcations. Theorem [T] states that the oriented subgraph 
Gv C G (illustrated in Figure [5J/bl) containing only vertex splittings, exists among 
gradient fields associated with convex bodies. Combinatorial expansion sequences 
([2]) associated with an n-vertex graph Q n 6 Q n appear on this oriented metagraph 
as an oriented path of length n — 2, starting at the root (Pz). Observe that in the 
kth step a vertex in the box-diagonal S + U = k is selected. Two such sequences are 
illustrated in Figure [B]/b2. Beyond theoretical interest, these metagraphs admit the 
study of interesting physical phenomena some of which we briefly discuss below. 

4.2. Inhomogeneous bodies. Our first example are inhomogeneous bodies. So 
far, throughout the paper we assumed convexity and material homogeneity. Re- 
laxing the latter constraint is equivalent to keep the convex surface bdK as the 
boundary of the body but let the mass be concentrated at the center of gravity C. 
As the location of C is varied in time as a curve vc(t), it generates a one-parameter 
family of gradient vector fields on bd K. A classical result in catastrophe theory 
states that the number of critical points of the gradient changes if and only if yq (t) 
transversely passes through one of the two caustics of the body [33]. Caustics 
(also known as focal surfaces) are the two surfaces formed by the curvature centers 
corresponding to the principal curvatures of bd K. Figure [Jj a)-c) shows the two 
caustics of an ellipsoid, a) corresponding to curvature minimum, b) correspond- 
ing to curvature maximum, c) intersection (superposition) of both caustics. When 
Tc(t) transversely crosses the caustic defined by the minimal principal curvature, a 
saddle and an unstable point meet at a saddle-node; when rc (t) transversely crosses 
the other caustic, a saddle and a stable point collide. Every saddle-node bifurcation 
corresponds to a vertex splitting (or face contraction, depending on the direction) 
on the quasi-dual Morse-Smale graph, so at each such event the corresponding path 
on G will move from one box-diagonal S + U = k to one of its neighbor diagonals. 
Figure [7J d) shows the different Morse-Smale graphs in the different domains de- 
termined by the intersections of the two caustics. It is easy to see that if C is 
located far enough from the center of gravity of the homogeneous body then the 
corresponding Morse-Smale complex is represented by the path graph P^. 

4.3. Collisional abrasion: chipping of rocks. Our second example is pebble 
abrasion via collisions. This process is most often described by averaged geometric 
PDEs, the most general such model is given by Bloore [7J as 

(5) v = 1 + 2bH + cK 

where v is the attrition speed along the inward surface normal, H, K are the Mean 
and Gaussian curvatures, b, c are constants. Solutions of these PDEs correspond 
to (inward propagating) wave fronts. The actual physical process is somewhat 
different: it is based on discrete collisons where small amounts of material are 




Figure 6. al) Metagraph Q corresponding to generic bifurca- 
tions in one-parameter families of vector fields on the sphere a2) 
Saddle-saddle connections inside the primary equilibrium class 
{S, U} = {2, 3} a3) Examples for convex bodies with some se- 
lected graphs bl) Oriented subgraph Q v c Q corresponding to ver- 
tex splittings b2) Two selected combinatorial expansion sequences 
b3) Morse-Smale complexes associated with real pebbles, derived 
with expansion sequences from the ellipsoid 
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being removed in a strongly localized area. Simple but natural interpretations of the 
discrete, physical abrasion process are chipping algorithms [17 , 34 and [29 where 
in each step a small amount of material is chipped off at point p by intersecting the 
body with a plane resulting as a small parallel translation of the tangent at p. We 
call such an operation a chipping event and their sequence a chipping sequence. 

In Section [3J we showed that any vertex splitting can be achieved by a suitably 
chosen convexity-preserving local truncation. It is not very difficult to show a re- 
lated, though converse statement: any sufficiently small chipping event will either 
leave the Morse- Smale complex invariant or result in one or two consecutive vertex 
splittings. If we regard the material abraded in a chipping event as a random vari- 
able x with very small, but finite expected value E{\) and very small variance (as it 
is often done in chipping algorithms) then most chipping events will be sufficiently 
small to form finite chipping sequences the subsequences of which are geometric 
expansion sequences of type ([3]). 

Chipping sequences do not represent a rigorous, algorithmic discretization of the 
PDEs, rather, they can be regarded as an alternative, discrete approximation of 
the phyisical process. As it was pointed out in [15] and [16], pebble surfaces display 
equilibria on two, well-separated scales. While the PDE description accounts for the 
evolution of global equilibria, local equilibria, corresponding to the fine structure 
of the surface are only captured by the chipping model. Figure [5] shows the high- 
accuracy scan of a real pebble with equilibria and the primal representation of the 
Morse-Smale complex. We can observe flocks of local equilibria accumulated around 
global equilibrium points. In fact, one motivation behind chipping algorithms is to 
better understand the interplay between the two scales. Chipping models appear 
to be successful in explaining laboratory experiments (cf. [29 j and 34 ) as well as 
geological field observations [3D ■ 

Figure|6]/b3 illustrates two, rather short geometrical expansion sequences leading 
to Morse-Smale complexes associated with real pebbles. While the abrasion of these 
pebbles was not monitored, given their simple Morse-Smale complex and nearly- 
ellipsoidal shape it is realistic to assume that these geometric expansion sequences 
are subsequences of the actual physical abrasion process (modeled by chipping 
sequences) which produced these shapes. Needless to say, many more experiments 
are needed to verify this theory. 



4.4. Concluding remarks. In this paper we showed that Morse functions A4c 
associated with convex bodies exhaust all combinatorial possibilities, i.e., to any 
two-colored quadrangulation Q* 6 Q of the sphere there exists a convex body K* 
such that Q* is the quasi-dual representation of the Morse-Smale complex associ- 
ated with K* . We proved our claim by showing that to each possible combinatorial 
expansion sequence Qi generated by vertex splittings, there exists a coupled geo- 
metric expansion sequence K% such that Qi is the quasi-dual representation of the 
Morse-Smale complex associated with Ki. Geometric expansion sequences were 
created by local, convexity preserving truncations of the convex body. Beyond 
proving the existence of all combinatorially possible convex shapes they also define 
a hierarchical order among them. We also showed that these sequences appear to 
be part of the natural geometrical description of collisional abrasion. 
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Figure 7. Caustics of an ellipsoid, a) Caustic denned by the 
minimal principal curvature b) Caustic denned by the maximal 
principal curvature c) Superposition of the two caustics d) Quasi- 
dual Morse-Smale graphs (in class Q) in the different domains, 
shown on the plane section through axes x and z 
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Figure 8. The Morse-Smale complex of a real pebble: primal 
Morse-Smale graph in class Q 3 * a) drawn on the pebble's surface 
b) drawn on the plane 



